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FREE-FLIGHT INVESTIGATION OF AERODYMMIC HEAT TRANSFER TO 
A SIMULATED GLIDE-ROCKET SHAPE AT MACH NUMBERS UP TO 10* 


By Andrew G. Swanson 

SUMMARY ^ 0<5'3d 

Heat-transfer measurements were made on a simulated glide-rocket 
shape in free flight at Mach numbers up to 10 and free-stream Reynolds 

numbers of 2 x 10^ based on distance along surface from apex and 3 x 10^ 
based on nominal leading-edge diameter. The model simulated the bottom 
of a 75° delta wing at 8° angle of attack. The data indicated that for 
the test conditions a modified three-dimensional stagnation-point theory 
will predict to reasonable engineering accuracy the heating on a highly 
swept wing leading edge, the heating being reduced by sweep by the 
3/2 power of the cosine of the sweep angle. The data also indicate that 
laminar heating rates over the windward surface of a highly swept flat 
glider wing at moderate angles of attack can be predicted with reasonable 
engineering accuracy by flat-plate theory using wedge local flow condi- 
tions and basing Reynolds numbers on lengths from the wing leading edge 
parallel to the surface center line. 


INTRODUCTION 


In view of the current interest in glide rockets, a free-flight 
investigation of aerodynamic heat-transfer characteristics of a glide- 
rocket shape has been made. The model configuration simulated a flat- 
bottomed delta wing having a blunt leading edge swept approximately 75° 
and flying at a trim angle of attack of 8°. In order to obtain a sym- 
metrical shape for flight test, the wing was simulated by a three-sided 
pyramidal shape, each surface of which would have heating characteristics 
similar to the undersurface of the glide missile; the edges of the pyra- 
mid would approximately simulate the leading edges of the glide missile 
wing. The test shape was tested with the same five-stage research mis- 
sile system used in the investigations reported in references 1 to 4. 

Although the model reached a maximum Mach number of 14.7 at an 
altitude of 88,100 feet, the stability of the fifth stage seemed to be 
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marginal and the model flew at a large and somewhat indeterminant angle 
of attack; thus, temperature data obtained during this portion of the 
flight are difficult to analyze. However, during the fourth-stage 
burning, heat-transfer data were obtained at Mach numbers up to 10 and 

free-stream Reynolds numbers of about 2 x 10^ based on the rearmost 

thermocouple location on the flat surface and 0.3 X 10^ based on leading- 
edge diameter. 

The fourth-stage rocket motor (JATO, 1.52-KS-33550, XM-19 (Recruit)) 
used in the present investigation was made available by the U. S. Air Force. 


SYMBOLS 


A 

A N 

Arp 

h 

k 

l 

M 

<1 

<*d 

R oo 

r d 


s 

t 


cross-sectional area at reference station, sq. ft 
normal acceleration, gravitational units 

transverse acceleration, gravitational units 
lift-curve slope per degree 

e nthalpy , Btu/ s lug 
thermal conductivity, Btu/ft-sec-9R 
distance between thermocouples, ft 
free-stream Mach number 
heating rate, Btu/sec-ft^ 
dynamic pressure, lb/sq. ft 

Reynolds number based on free-stream conditions 

Reynolds number based on free-stream conditions and twice 
leading-edge radius 

surface area used in conduction corrections, sq ft 
time, sec 

temperature, °R unless otherwise noted 


T 
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w 

weight, lb 


X 

distance, ft 


dU 

rjy 

velocity gradient, l/sec 


a 

angle of attack, deg 


e 

orientation of resultant acceleration vector, deg 


A 

sweep angle, deg 


P 

viscosity, slugs/ft-sec 


P 

density, slugs/cu ft 


Subscripts : 


theo 

theory 

/- . 

w 

wall 


s 

stagnation 


3D 

three-dimensional 


aw 

adiabatic wall 


00 

free stream 


res 

resultant 


n 

thermocouple number 



MODEL AND TEST 



Model Configuration 

.•. s 

The fifth stage of the missile system used in this investigation 
consisted of a three-sided pyramidal test nose which had a blunted apex, 
a stepped cylindrical midsection, and a 20° total angle conical frustum 
tail. A photograph of the fifth stage is shown in figure 1; pertinent 
dimensions are given in the sketch (fig. 2) which also shows dimensions 
of the forth- and fifth-stage combination. The test nose simulated a 
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glide rocket shape, and nose details are shown in the sketch of figure 3 
and the photographs of figure 4 . 

The test nose was fabricated from Inconel of 0.05-inch nominal 
thickness. One edge of the pyramid was rolled to an exterior radius of 
0.10 inch, and the resulting two surfaces were welded at the edges to a 
third surface. The pyramid thus formed was welded to an Inconel cylin- 
der at the lines of intersection of the pyramid and cylinder. The for- 
ward end of the nose was formed from a block of nickel. 

A conical radiation shield was mounted inside the test nose and 
bolted to the nickel block at the forward end. The after end of the 
test nose was slip-fitted over a micarta block attached to the radiation 
shield. The thermocouple switching motor was attached to a bracket 
within the radiation shield. 

A double-walled Inconel cylindrical section behind the nose housed 
most of the telemetering equipment. The stepped conical section con- 
tained an alumina insulating ring, which was used to isolate electri- 
cally the two ends of the fifth stage to form the telemeter antenna. 

The Fiberglas lined, Inconel conical frustum at the after end was intended 
to give the fifth stage static stability and also to serve as an extension 
to the rocket-motor nozzle . 

The pyramidal nose was polished to a general surface roughness of 
10 microinches or less as determined optically by a fringe-interferometer 
microscope. There were, however, several pits and scratches which were 
considerably deeper than the 10-microinch level; these did not seem to 
have significant effects on the data measured. 



Instrumentation 

Measurements of the output of 24 thermocouples and 4 accelerometers 
were transmitted from the test vehicle during flight by a six-channel 
telemeter. The chromel-alumel thermocouples were made from Wo. 30 gage 
wire and were spot welded to the inner surface of the skin of the test 
nose . The two leads of each thermocouple were welded separately to the 
skin with a spacing between the leads of about 1/32 inch to l/l6 inch. 
Those thermocouples located behind the leading-edge stagnation line had 
the spacing lengthwise along the leading edge; that is, both thermo- 
couple leads were on a line parallel to the stagnation line . The thermo- 
couple located in the nickel block at the forward stagnation point was 
welded at the bottom of a tapped hole and the leads were led through the 
hollow bolt which supported the radiation shield. The locations of the 
thermocouples are shown in figure 5- Wall thicknesses measured at the 
various thermocouple stations are given in table I. 
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Three constant voltages and outputs of 12 thermocouples were com- 
mutated by a switching motor and transmitted by one telemeter channel 
at a sampling rate of about 6 per second; the remaining 12 thermocouples 
and the 3 constant voltages were transmitted in similar fashion on 
another telemeter channel. The constant voltages served as in-flight 
references for calibration of the thermocouple-telemeter system. The 
voltages were chosen to be equivalent to the low end, middle, and high 
end of the range for which the thermocouples were calibrated. Since 
the range was chosen on the basis of possible turbulent heat transfer 
and since the flow was primarily laminar during the test, the actual 
measured temperatures were less than one -half of the full-scale range . 

Each of the four remaining telemeter channels was used to transmit 
a continuous measurement of the output of an accelerometer; two of these 
measured longitudinal accelerations, one measured normal accelerations, 
and one transverse accelerations . The accelerometers were calibrated in 
standard earth gravitational units g for the following ranges: 

Longitudinal accelerometers -50g (drag) to 150g (thrust) 

Normal and transverse 

accelerometers -25g to 25g 

Other instrumentation consisted of ground-based radars for measuring 
velocity during the early portion of the flight and for determining space 
position of the missile system. A radar- tracked radiosonde balloon was 
used to determine atmospheric conditions and wind velocity; these data 
were determined at the altitude of the high-speed portion of the flight 
within about one-half hour of the flight . 


Flight 

The propulsion system consisted of five stages of solid fuel rocket 
motors. The test nose, which was attached to the fifth stage, and the 
booster stages are shown on the launcher in figures 6 and 7- Character- 
istics of the rocket motors and stage weights for a similar five-stage 
system are tabulated in reference 3; these data are also applicable to 
the test reported herein. 

The missile system was launched at an elevation angle of 73° • The 
first two stages were used to propel the remaining three stages to a 
peak altitude of about 9^->200 feet. Just after apogee, when the flight 
path was inclined downward at about 2°, a preset mechanical timer fired 
the third stage. Shortly after third-stage burnout, a delay squib, 
ignited at third-stage ignition, fired the fourth stage; firing of the 
fifth stage was accomplished by use of a pressure switch, mounted on 
the fourth stage, which closed when the chamber pressure decreased as the 
fourth-stage motor burned out. 
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DATA REDUCTION 


Trajectory Data 

Velocity of the missile system for the first 28 seconds of flight 
was obtained from the ground-based Doppler radar. Thereafter, velocity 
was obtained from differentiation of the space-position radar data until 
third-stage ignition; from then until fifth-stage burnout, velocity data 
were obtained from integration of the longitudinal accelerometer data. 

The flight path was determined from the space-position radar data from 
launching time until third-stage ignition; shortly after this time the 
radar ceased to track the missile system. The flight path was then recon- 
structed by double integration of the longitudinal accelerometer data. 


Heat-Transfer Data 

The aerodynamic heat transfer to the model was obtained from the 
thermocouple data which gave temperature time histories for the inside 
surface of the skin of the test nose. These temperature data are nor- 
mally reduced by use of the reference calibration voltage measured after 
each sampling cycle of the thermocouple. In this test, shifts were noted 
in the reference voltages at the times of rocket-motor firing and small, 
slow, but persistent drifts at other times. Consequently, the thermo- 
couple data were reduced by use of calibration voltages recorded both 
before and after each sampling cycle of the thermocouple outputs . These 
points were plotted and a curve drawn through the average of, or the 
most consistent of, the data points obtained at each time. 

These faired inside-surface temperature curves, together with the 
one-dimensional heat-flow-analysis method of reference 5 for thermally 
thick walls, were used to compute outside wall temperatures . This com- 
puting method exhibits tendencies toward instability when small time 
intervals are used if the time rate of increase of inside wall tempera- 
ture is small. In order to obtain more points at a larger time interval, 
two separate computations of outside temperature were made at intervals 
half as large and then the two were joined together to give the outside 
temperature history. The analysis assumes constant thermal properties 
for the wall, and the values used were those for a temperature midway 
between the highest and lowest measured temperatures for each station; 
this use of constant thermal properties is believed to introduce small 
error in the final heat-transfer data . In general, the computed out- 
side wall temperatures were used directly in the one -dimensional analysis 
method of reference 5 "to determine the aerodynamic heat input to the skin 
of the test nose, again on the assumption of constant wall thermal 
properties . 

■s 
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For the thermocouples on the leading edge, the actual physical 
thickness of the skin was not used in computing temperature and heating 
rates. Since the heat flow into the outer surface tends to be "focused" 
toward the center of the leading-edge radius, an effective leading-edge 
thickness equal to the ratio of the leading-edge volume to the external 
surface area was used. For the geometry of the configuration tested, 
the effective thickness was about three-quarters of the actual thickness 
The heating rates computed with this effective thickness were about 
10 percent to 30 percent less than the rates computed with the actual 
leading-edge thickness . 


Corrections for Conduction 

Conduction effects in the longitudinal direction (or along any ray 
from the apex) on the test nose were estimated to be negligible . Con- 
duction effects in the lateral direction (or across the nose) were also 
estimated to be negligible except near the corner or leading edge . Lat- 
eral conduction corrections were approximated by a method similar to 
that of reference 3 • A lateral strip of unit width on the nose was 
divided into blocks approximately centered about a thermocouple . The 
change in heat flow at a station n is then given by 

An = ^n+l T n - T n+1 , ^-1 T n ~ T n-1 
n ^n ^n+1 ®n ^n-1 

At the corner, the block surface area S n was taken as the area 

between points where the leading-edge circumference became tangent to 
the flat surfaces; the remaining block edges were then taken at stations 
midway between thermocouples. The areas A n are, of course, the cross- 

sectional area of the skin; lengths between thermocouples Z n were 

taken at a distance two-thirds of the way out from the inside surface . 
Outside surface temperatures were used in these calculations . 

This method of correcting for lateral heat flow is not exact; it 
assumes, for example, that the heat flows through and along the skin 
can be calculated independently and the resulting solutions superimposed 
It is believed that insufficient temperature-distribution data were 
obtained to warrant a more rigorous analysis. The general conclusions 
regarding the overall heat transfer would not, however, be likely to be 
modified by even fairly large changes in the estimate of these lateral 
heat flow corrections . 

Radiation heat losses were estimated to be negligible and were 
therefore not included in the analysis. 


• • • • • 
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ACCURACY 


Trajectory Data 

The possible inaccuracies in the Mach number at the time of third- 
stage ignition are estimated to be less than ±0.1. Integration of the 
accelerometers for velocity after this time gives an accumulative error 
with time. This would result in an additional possible error in the 
Mach number at fifth-stage burnout which is estimated to be less than ±0.5- 
In view of the consistency between preflight trajectory estimates and the 
flight test results for velocity increments for each stage of the flight 
for this test and the tests reported in references 1 to if (for which simi- 
lar propulsion systems were used), the Mach numbers are believed to be 
more accurate than the foregoing figures would indicate. 

Uncertainties in the flight-path angle due to failure of the space- 
position radar to track the test vehicle after third-stage ignition 
result in estimated possible inaccuracies in altitude at fifth-stage 
burnout of less than ±2,000 feet. 


Temperature and Heat-Transfer Data 

The full-scale temperature range (chosen on the basis of possible 
turbulent flow at M = 15) was 1700° F, which results in an estimated 
possible basic temperature data inaccuracy of ±17° F. The level of the 
curve faired through the measured data points has an estimated possible 
error of about ±5° F. Since the heating rates were fairly low during 
the primary-data-measurement period (fourth-stage burning) possible 
inaccuracies computed for the final heat-transfer data would be fairly 
large . However, in view of the reasonable consistency of most of the 
data, the accuracy is believed sufficient to warrant the conclusions 
drawn. A more quantitative picture of the accuracy of the heat-transfer 
data is presented in the "Results and Discussion" section. 


RESULTS AND DISCUSSION 


Flight Path 

Flight-test data were obtained as the test vehicle moved along the 
trajectory shown in figure 8. Times of significant events are indicated 
on the trajectory. Atmospheric condition^ determined from the radiosonde 
data are shown in figure 9 as functions of time from third-stage ignition 
until fifth-stage burnout. Altitude and velocity data are shown in fig- 
ure 10 for the same time interval. The telemeter signal from the test 


NACA RM L58G03 


9 


• • ••• • • 



% 


vehicle failed momentarily just before fifth-stage burnout and ceased 
abruptly just as the fifth stage began to decelerate after burnout; 
evidently the test vehicle failed structurally at this time . 


Test-Vehicle Stability 


Plotted in figure 11 are the magnitudes of the normal and transverse 
accelerations measured by the accelerometers, and in figure 12 are the 
angles of attack estimated from these accelerations . The accelerations 
and estimated angle of attack were of negligible magnitude up to the time 
of fourth-stage ignition. The fourth- and fifth-stage combination evi- 
dently received a disturbing force as it separated from the third stage . 
This disturbance was fairly well damped near the time of fourth- stage 
burnout . 

At fifth-stage ignition (t = 89.5 sec), the test vehicle started a 
somewhat divergent motion; the motion was not, however, purely divergent, 
as is shown by the oscillatory nature of the acceleration data (fig. 11) . 
The exact reason for this apparent instability is not known; probably 
there was insufficient static stability, although it is possible that 
the vehicle was damaged at separation from the fourth stage. The test 
vehicle was originally developed for tests of blunt noses; the calculated 
stability for the sharper (and higher lift) nose of the present test 
was less than for the blunter nose models, although still presumably 
adequate. Probably the stabilizing force of the flare was overestimated 
and the use of the sharper nose changed the amount of stability from just 
adequate to rather marginal . Note that tests of the blunt-nose configu- 
rations of references 1, 3, and 4 showed small or negligible angle of 
attack while the test of a sharper nose of reference 2 showed significant 
angle of attack (although less than the angles apparently reached by the 
test vehicle of this test) . 

As is shown in the angle- of -attack data (fig. 12), there was a 
tendency for the fifth stage to trim at angle of attack; the steadily 
decreasing weight as the propellant burned would result in the measured 
increasing normal acceleration for a constant trim angle. There was 
undoubtedly some rolling motion; however, polar plots of the resultant 
acceleration vector established no consistent trend to a rolling motion, 
and it is believed that the motion was mainly purely oscillatory in the 
pitch and yaw planes. 


From the resultant acceleration vector, 
were computed from the equation 


a. 


res 


w/a n 2 + At 2 
c Iu qdA 


resultant angles of attack 
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where linear weight variations with time were assumed during rocket- 
motor burning. Since force data were not available for the configura- 
tion tested, estimates of zero-lift C]^ were used, 0.12 per degree 

based on a 9-inch-diameter body area for the fourth- and fifth-stage 
combination and 0.07 per degree based on the 6.2-inch-diameter body 
cross-section area for the fifth stage. This probably gives an over- 
estimate of the actual model angle at the higher angles of attack, since 
the lift-curve slope would be expected to be nonlinear and to increase 
with angle of attack. 

The angles of attack plotted in figure 12 are the angles determined 
from the resultant acceleration vectors and are average angles faired 
through the oscillations. Also plotted in figure 12 are the orientations 
of the resultant acceleration vector. The motion can be seen to be one 
primarily in the pitch plane. Note that only in the time from about 88.4 
to 89.5 seconds (the latter part . of fourth-stage burning) are the esti- 
mated angles of pitch and yaw small (l° or 2° or less) . Contributing 
to the large angle of attack occurring at fourth-stage ignition is the 
relatively low dynamic pressure at this time (about 500 pounds per square 
foot); this dynamic pressure steadily increased to about 2,700 pounds per 
square foot at fourth-stage burnout and about 6,200 pounds per square 
foot at fifth-stage burnout. 

The approximate angles of attack calculated are referred to in the 
rest of the discussion as estimated angles; the limitations in accuracy 
discussed above are implied. 


Basic Temperature Data 

The values of inside surface temperatures determined from the telem- 
eter records are shown for several typical stations in figure 13 for 
times from 87 seconds on. The data before t = 87 seconds had similar 
scatter. Shown in figure 13(a) are data for thermocouples 2 and l6 
(located on the leading-edge stagnation line) and in figure 13(b) are 
data for thermocouples 4 and 21 (located on the center line of the flat 
surface which contained most of the thermocouples) . The lower set of 
curves on each of these figures was obtained from one thermocouple chan- 
nel; the upper set of curves was obtained from the other thermocouple 
channel. Faired through the data points (which are shown as reduced from 
calibrate voltages before and after each sampling cycle as mentioned in 
the data reduction section) are the curves used for computation of out- 
side surface temperature. Shown in figure 13(c) are the data points and 
faired inside surface temperature for thermocouple 1, which was located 
behind the nose stagnation point, along with a computed outside wall 
temperature variation. Thermocouple 1 was installed to give a qualita- 
tive indication of possible nose-tip melting; uncertainties as to 
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conduction corrections and directions of heat flow into this thermocouple 
precluded any analysis of heat transfer at this station. It is apparent, 
however, that high heating rates existed at the stagnation point and that 
at about the time of fifth- stage burnout, the forward stagnation surface 
was beginning to melt; this melting may have contributed to the model 
failure which evidently occurred as the test vehicle began to decelerate. 
The temperatures were of the order of those that would be exnected for 

the flat-faced nose stagnation point. 

* 

Inside temperatures read from the curves faired through the measured 
data points are tabulated for all thermocouples and for several times 
during the flight in table II. 


Computed Outfcide Wall Temperatures and Heating Rates 

Shown in figure 14 are typical time histories of computed outside 
wall temperatures . The data are presented for the time periods during 
which the angle of attack was sufficiently small and the heating rates 
sufficiently high to provide reasonable heat-transfer data. Shown are 
the faired inside temperatures, the computed outside temperature points, 
and the curve faired through these computed points. The outside tem- 
peratures computed could generally be used directly in the computation 
of heat-transfer rates; however, stable oscillations arising in the 
numerical analysis occasionally necessitated fairing. The heating rates 
were in some cases sufficiently high for significant temperature gradi- 
ents through the skin to exist. 

The computed heating rates to the outside surface (one-dimensional 
heat flow through skin with lateral conduction neglected) are shown for 
several typical stations in figure 15 . Both computed points and the 
curve faired through them are shown. When the actual aerodynamic heating 
rates are low, slopes of the measured temperature time histories, which 
are proportional to heating rate, are, of course, also low. At these 
times inaccuracies in temperature slope determination result in fairly 
large amounts of scatter in the heat-transfer data. 

Also shown in figure 15 are typical data including conduction cor- 
rections . The corrections are by no means negligible for stations on 
or near the leading edge; however, fairly large percentage changes in 
the conduction corrections would be required to alter the basic con- 
clusions drawn from the data. 

An idea of the possible quantitative accuracy of the final heat- 
transfer data is probably best determined ' from the scatter of the data 
of figure 15 - 
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Theory Used for Leading-Edge Heat Transfer 

The theory used for comparison with the experimental leading-edge 
heat transfer was a modification of the Fay and Riddell theory (ref. 6) 
for three-dimensional stagnation-point heating. This theory was used 
mainly because of its ease of calculation (when appropriate simplifica- 
tions are used) . The basic theory of reference 6 is 

^theo^ = °-^(P w ^w) 0 ' 1 (Ps^s) 0 ' 1+ ( h s “ h w)(d^) s> 3D 


for a Prandtl number of 0.71 and a Lewis number of 1 . 0 . Reference 7 
shows that, for the conditions of this test, varying Lewis number 
between 1 and 2 has small effect on the heating rate. The above expres- 
sion is quite easily evaluated if the Sutherland variation of viscosity 
with temperature is used (ref. 6 indicates this law to be reasonably 
valid for temperatures up to 9000° K) and if ideal gas relations are 
used. The use of ideal-gas rather than real-gas relations should result 
in negligible differences in heating rates for the conditions of this 

test (as indicated in ref. l) except possibly for the effect in ( ^ 

\dx 

This parameter was evaluated from the theory of Korobkin for Mach numbers 
up to 5 and from Newtonian flow pressure distributions for higher Mach 
numbers. (See curve in fig. 16 of ref. 1.) The use of real-gas rather 

than perfect-gas relations for ( — ) would lower the theoretical heating 

\dx/s 

rates by about 5 to 10 percent at M = 10 and by lesser amounts at lower 
Mach numbers . It can be argued that for application of the theory to the 

swept leading edge, the real-gas effects on ( — ) should be applied on 

\dx/s 



the basis of Mach number normal to the leading-edge shock rather than 
free-stream Mach number; for the conditions of this test, there would 

then be negligible real-gas effects on f— ) • 


Two additional factors were necessary to convert this three- 
dimensional heating rate to a rate for a swept leading edge. First, 
the three-dimensional rates were converted to a two-dimensional rate 

by multiplying by the factor which was used by Lees (ref. 8) and 

r2 

was shown to be approximately correct by Reshotko and Cohen (ref. 9 ) • 

T 

For the ratio — & for this test, reference 9 would indicate that this 

T s 

factor might be more appropriately 0.72 to 0.73 instead of O. 707 . Sec- 
ond, this heating rate at the stagnation point of an unswept cylinder 
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must be modified to account for the reduction in heating rates due to 

sweep. For the results given herein, the factor used was (cos A)-'*/^; 
the validity of this factor (which has some theoretical justification) 
is discussed subsequently. 


Experimental Leading-Edge Heat Transfer 

The experimental heating rates on the leading (or pyramid) edge are 
shown in figure 16 as ratios with respect to theoretical heating rates 
for several times (and Mach numbers) during the flight. Data are pre- 
sented for times near the end of third-stage burning when the measured 
rates were somewhat greater than the scatter in the data and for times 
during f earth -stage burning when the estimated angle of attack was rea- 
sonably low. 

The experimental data, in general, are in agreement with or are 
less than the values predicted by the theory given in the previous sec- 
tion. The data for 86 seconds and 87 seconds were obtained during third 
stage burning when there was marked scatter in the computed heating rate 
points (fig. 15)- The agreement with theory may be fortuitous under 
these conditions; however, the four measuring stations on the leading 
edge do show remarkable consistency. The reason for the low ratios 
shown for the data at 88. A through 89. 0 seconds cannot be ascertained 
with certainty. Possibly, the angles of attack and yaw are still suf- 
ficiently large to affect the measurements . As is shown subsequently, 
the angle of attack, although estimated to be fairly low, does have a 
noticeable effect on the relative levels of heating on the three flat 
faces up to about 89 seconds. The use of estimated angle of attack 
(modifying effective sweep angle) in computing theoretical values of 
leading-edge heating rates, however, had negligible effect at times 
after 88.4 seconds. 

There are probably Mach number effects in the reduction with sweep 
parameter (the theory of ref. 10, for example), which are not accounted 
for by the method of prediction used herein. These effects would, how- 
ever, be expected to raise the predicted level at the lower Mach numbers 
or increase the difference between prediction and experiment. 

The data for times after 89 seconds show quite good agreement with 
the theory. At these times the measured heating rates are highest and 
the data would be expected to have the most accuracy. 

The data of figure 16 would indicate that the reduction in heat 

transfer with sweep is proportional to (cos A)^/^. If cos A (which 
has also been proposed as the proper value for this parameter) had been 
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used instead of (cos Ap/ 2 in computing the theoretical heating rates, 
the level predicted by the theory would have been 1.93 times higher than 
the theory level shown in figure 1 6. For the data of this test, the use 
of cos A would therefore result in a large overestimate of the leading- 
edge heating rate, (it should be noted that the difference between 

cos A and (cos Ap/ 2 is of appreciable magnitude only at the higher 
values of sweep angle, say greater than 40°.) 

From the data of the test, it would appear that the heating of a 
swept cylindrical leading edge of a glide missile wing can be predicted 
with reasonable accuracy by using the Fay and Riddell theory as modified 
herein. Since this theory is very easy and straightforward to evaluate 
(when Lewis number is assumed to be 1.0 and when perfect-gas relations 
are used), it would seem to be a useful relation for engineering calcu- 
lation of leading-edge heating. The discussion presented in the fol- 
lowing section should, however, be considered in an evaluation of its 
merits . 


Factors Arising from Model Geometry 

Before comparing the data of this experiment with other experiments 
and theories, it is desirable first to consider some factors peculiar 
to this test. The leading edge of the model is not a true cylindrical 
leading edge in that it is formed by the radius at the intersection of 
two sides of the pyramidal nose. (See fig. 3-) The precise effect of 
this geometry on the leading-edge stagnation-line heat-transfer rates 
is difficult to evaluate. Since the sonic line on a cylinder is about 45° 
from the stagnation line, and since the point of tangency of the side 
and leading edge of the configuration tested was 60° from the stagnation 
line, it would appear that the lack of a complete cylinder at the leading 
edge would not significantly alter the test results from those that would 
be expected for a swept infinite cylinder. However, consideration of 
the shock structure around the test nose complicates the analysis . 

This shock structure is difficult to determine without visual flow 
tests of the configuration flown. A rough qualitative analysis can be 
made, however. The shock over the 8° flat surfaces of the simulated 
wing would, of course, be detached from the surfaces. The juncture of 
these shocks at the edges of the surfaces (the simulated wing leading 
edge) would probably alter the shock around the leading edge and give 
it somewhat different curvature and detachment distance than would 
exist on a yawed infinite cylinder. Since the flow-deflection angle 
is low, it is believed that this effect would be small, particularly 
at the higher Mach numbers when the shock would lie quite close to the 
flat surface . The shock formation over the nose can also be considered 
from another point of view. It can be reasoned that the shocks over 


the simulated leading edges would he approximately the same as the coni- 
cal shock occurring on a cone circumscribed about the pyramid, or a 
13-5° half-angle cone. The shock over this cone would also lie close 
to the surface, particularly at the higher Mach numbers, and its effects 
on changing the heat transfer from the value that would be expected on 
a yawed infinite cylinder are believed to be small. 

The effect of thse shocks would probably be to cause variation in 
heat transfer along the length of the leading edge. (The blunt-nose tip 
would induce a normal shock at the nose which would also tend to produce 
the same effect.) The data of figure 16 show small, and no consistent, 
variations with length along the leading edge; the small differences in 
heating rates between the measuring stations at any given time are within 
the overall accuracy of the data. This fact, plus the generally good 
agreement with the theory, is believed to warrant the above statement 
that the shock formation over the test nose does not significantly alter 
the heating rates measured from those that would occur on a swept infi- 
nite cylinder. The argument that the agreement between theory and experi- 
ment may be fortuitous cannot, however, be rigorously disproved. 

The use of an effective leading-edge thickness less than the physi- 
cal thickness in the data reduction for the leading edge (see data reduc- 
tion section) reduced the computed heating rates from those computed by 
use of the actual physical thickness. If the actual thickness were used, 
the experimental rates would be about 30 percent higher at the higher 
values of time presented and higher by somewhat lesser amounts at the 
earlier times; the change would be negligible for the data presented 
before 88. k seconds. The use of the effective thickness rather than the 
actual thickness is considered the more reasonable procedure. The agree- 
ment between theory and experiment would still be fairly reasonable, how- 
ever, if the thicker wall were used in the data reduction. 


Comparison With Other Theories and Experiments 

Theroetical heating rates were also computed by the theory of ref- 
erence 11 (where the parameter — (eq. 10(a), ref. ll) was chosen 

\ \Tc 

to be 1.0 which is equivalent to assuming that conductivity and viscos- 
ity of air vary in the same way with temperature ) . Heating rates were 

computed by using both =1.0 and — i — = O.Sj. (This value of 

T s T s 

"recovery factor" for a sweep angle of 75° was obtained from ref. 12.) 
The heating rates thus computed were, respectively, slightly higher than 
and slightly lower than the rates computed from the modified Fay and 
Riddell theory. The theory of reference 11 predicts, for M = oo (or 
practically for M 7) a variation of heat-transfer coefficient with 
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sweep proportional to (cos A) ; the theory predicts a lesser reduction 
with sweep for lower Mach numbers (as does the theory of ref. 10) . The 
experimental results of reference 11 also show a reduction in heat-transfer 

3 /2 

rates due to sweep proportional to (cos A) ' ; the data were obtained, 
however, at a maximum sweep angle of about 45° • 

The experimental data of reference 13 also indicate that the heating 

rate is decreased by (cos A)^/ 2 for sweep angles up to 45°. (There is 
considerable scatter in the data, but the general trend is a reduction 

proportional to (cos A)V 2 .) However, the data for a 70° sweep angle 
show a reduction more nearly proportional to cos A. The reason for this 
change in trend is not readily apparent. It should be noted that the 
data of reference 13 were obtained at low free-stream Reynolds numbers 
(315 to 2,100) . 

Experimental data reported in reference 14 indicate that the reduc- 
tion in heating rate with yaw angle is proportional to (cos A)^'^. The 
data were obtained at a nominal Mach number of 10.4 and at sweep angles 
up to 70°- However, these test results are subject to a possible limita- 
tion in that they were obtained at Reynolds numbers so low that the 
possibility of the occurrence of slip flow during the test should be 
considered. 

The experimental data of reference 12 indicate a reduction with 
sweep in heat-transfer coefficient which is proportional to cos A. The 
reason for the reduced benefits of sweep in reducing leading-edge heating 
which is indicated by these data is not known. At least part of the dis- 
agreement is attributable to the fact that the data of reference 12 are 
heat-transfer coefficients. Heating rates would show a greater benefit 
of sweep in alleviation of heating than do heat- transfer coefficients 
because of the influence of recovery factor variation with sweep. Recov- 
ery factor decreases with sweep, hence, T aw and the resulting forcing 

function T aw - T w also decrease with increasing sweep angle. 

The theory of reference 10 also shows that, for high Mach numbers 
(M = 7), sweep reduces the heat- transfer coefficients by a factor of 

about (cos A)^/^. As was previously mentioned, this theory (and others) 
also shows a Mach number effect on the reduction with sweep parameter; at 
lower Mach numbers the effect of sweep is theoretically not so beneficial. 
The experimental data of the test reported herein show different trends. 

The reason for this difference is not readily apparent; it might possibly 
lie in inaccuracies in the data at the lower Mach numbers (M « 4 at 86 
and 87 seconds) and incorrect estimates of angle of attack for the data 
from t = 88.4 to 89. 0 seconds (fig. l6) . 
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Comparison With Theory for Data on the Flat Wing Surface 

The experimental heating data on the flat surface of the glide mis- 
sile wing are compared with theory in figure 17 . The data are shown as 
reduced to local Stanton number. Local flow conditions were those com- 
puted for an 8° wedge and by use of total pressure losses through an 
oblique shock at the wedge leading edge . Real-gas effects on total tem- 
perature and specific heat were used, the viscosity was assumed to be 

proportional to t ^-76 an ^ ideal-gas relations were used to compute 
the remaining flow conditions . 

Theoretical Stanton numbers were computed from the theories of 
Van Driest (ref. 15 ). Adiabatic wall temperatures were computed by 
using recovery factors of 0.84 and 0.88 for laminar and turbulent flow, 
respectively. The laminar recovery factor was used in the reduction of 
experimental data. 

The experimental data of figure 17 are shown as correlated by two 
different methods. In method 1 (figs. 17(a) to (c)) the data are plotted 
on the basis of length from the leading edge to the measuring station in 
a direction parallel to the wing center line. In method 2 (figs. 17(d) 
and (e)) the data are correlated on the basis of length along the rays 
from the model apex. There is, in general, better correlation by the 
first method (figs. 17(a) to (c)). The data indicate that the flow over 
the flat surface was probably laminar. 

Data are shown in figure 17 for measuring stations on all three 
flat surfaces. The data for the sides opposite the main thermocouple 
array, shown as flagged symbols, scatter about the general level of the 
other data. The single flagged symbols are for thermocouples 13 and 23 
and the double flagged symbols are for thermocouples 14 and 24 (see 
fig. 5 ) . Time history comparisons of the data for thermocouples on all 
three sides are shown subsequently. 

There is scatter in the data of figure 17, the points on or nearest 
the center line showing the greatest variations . Although it is prob- 
able that the simplified flow analysis used does not properly account 
for the effects of the actual pressure distribution on the heat transfer, 
it is possible that the anomalies in the data are due (at least in part) 
to inaccuracies in the basic data. The points showing the greatest 
scatter (both from the general trend at any one time and in self- 
consistency at the various times) are, in general, at measuring stations 
having the lowest heating rates . The exception to this trend is the 
point farthest forward and outboard on the surface (thermocouple 3) • 

The problem at this station may be one of conduction correction. Since 
adequate temperature distributions were not obtained at the longitudinal 
station on which this thermocouple was located, the same correction was 
applied to thermocouple 3 that was applied to thermocouple 9 • 
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The points showing the most deviation from theory are at t = 88.4 
and 88.6 seconds, where there is still some angle of attack. The use of 
the estimated angles of attack in computing local flow conditions, how- 
ever, shows small effect on the heat transfer due to angle of attack 
(as was true for the leading-edge data shown in the previous section) . 

It does appear, however, at least for the conditions of this test, 
that the heat transfer to the windward surface of a highly swept glide 
rocket wing at moderate angles of attack can he predicted with fair 
accuracy by the use of the Van Driest theory with Reynolds number based 
on distance from the leading edge parallel to the wing center line 
(method 1), and using appropriate local flow conditions. 

It is interesting to note that, for this test, the computed laminar 
heat- transfer rates on the wing surface were not much different (about 
5 percent less) when local flow conditions were based on losses through 
a normal shock at the wedge leading edge instead of assuming an oblique 
shock at the wedge leading edge. Therefore, for this test the agreement 
between theory and experiment would not change significantly if either 
type of local flow conditions were used. 

At t = 91 seconds (or M = 14.4), laminar Stanton numbers were 
calculated for the stations on the windward surface by the method of 
Van Driest (ref. 15). With local flow conditions computed as discussed 
previously (but with the wedge angle taken to be 8° plus the angle of 
attack of approximately 12°), theoretical heating rates were computed. 
These theoretical rates were of the same general magnitude as the meas- 
ured rates. This agreement may very well have been fortuitous, since 
there is some uncertainty as to the validity of the angle-of-attack 
determination, and these data have therefore not been plotted. 


Comparison of Data on the Three Faces 

In addition to the main array of thermocouples on one face, addi- 
tional thermocouples were placed at two longitudinal stations on the 
center line of each of the remaining two surfaces . The data from each 
of the three surfaces are compared in figure 18, both outside wall tem- 
perature and heat-transfer rates are shown for the portion of the flight 
commencing with third-stage ignition. 

Up to the time of fourth-stage ignition (t = 87-4 sec) the data are 
in close agreement as would be expected; differences are due mainly to 
scatter in the basic data due to the low temperatures and heating rates 
during this time. 

The angle of attack which occurs at fourth-stage ignition 
(t = 87A sec) affects the heating rates in a logical manner, the 
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windward side (measured by thermocouple Ik) of the model showing the 
highest rates . There are noticeable differences between the windward 
side and the two leeward sides (measured by thermocouples 13 and 15 ) . 

At about 88.8 to 89 . 0 seconds the heating rate on the windward side begins 
to decrease (as the angle of attack is small at these times) and approaches 
the level of the data for the remaining two sides. At about 89.5 seconds 
at fifth-stage ignition when the test vehicle again goes to a high angle 
of attack, the data for the three sides of the model again diverge and 
again the trends are as would be expected. 

When the data of figure 18 are considered in comparison with the 
data of figure 12, note that the model has accelerated from a Mach num- 
ber of 6 to a Mach number of 10 during the time interval from maximum 
angle of attack (t 555 88 sec) to the time of maximum heating on the wind- 
ward face (t « 89 sec) . The maximum heating rates, therefore, would 
not be expected to occur at the same time as the maximum angle of attack. 


CONCLUDING REMARKS 


A free-flight rocket-model test of a simulated glide rocket indi- 
cates that for Mach numbers up to 10 and for the free-stream Reynolds 

number range of the test (2 X 10^ at the rearmost thermocouple location 
on the flat surface and 0.3 x 10 based on leading-edge diameter): 

1. Stagnation-line heat-transfer rates to a highly swept cylindrical 
leading edge of a glide missile wing can be approximated with reasonable 
engineering accuracy by the use of the three-dimensional Fay and Riddell 

theory modified to two-dimensional heating rates by the factor — and 

reduced by the three-halves power of the cosine of the sweep angle. 

(This latter factor is not in agreement with reductions suggested by 
some wind-tunnel results, which have indicated that the reduction in 
heating with sweep varies as the cosine of the sweep angle.) 

2. For the highly swept flat surfaces of the glide vehicle investi- 
gated, the laminar heat transfer to the windward surface at moderate 
angles of attack can be approximated with reasonable engineering accuracy 
by the use of the laminar flat-plate theory of Van Driest by using wedge 
local flow conditions and basing Reynolds numbers on lengths from the 
surface leading edge parallel to the surface center line (rather than on 
length along rays from the apex of the swept surface . 

Langley Aeronautical Laboratory, 

National Advisory Committee for Aeronautics, 

Langley Field, Va., June 27, 1958* 
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Figure 1.- Photograph of test nose mounted on fifth stage. 
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(b) Missile system mounted on launcher just prior to elevation for firing .L-58-I698 


Figure k . - Photograph of test nose . • 
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Figure 5.- Sketch of thermocouple locations on nose. Dimensions in inches. 
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Figure 6.- Missile system being prepared for launching. 
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Figure J.- Missile 


system elevated on launcher ^“57-3841 


CONFIDENTIAL 


100x10 


MCA RM L58G03 


31 


• • ♦ • • • • • 



• • 




• • 


• • • • • 

• • • • ••• • 



-P 

cti 

Pn 

-P 

b 0 
•H 
rH 
<H 

0) 

I 1 

0 

■H 

-C 

•p 

OJ 

0) 

EH 

1 

co 


b0 

•H 

PH 


* •pni'HtY 



icao 


32 


• • • • 
• • 

• • • 
• • 

• • • • 


• • • • 


• • • • • , 


MCA RM L58G03 



Si. 


s 

•H 

E-< 


& 


09B/4J 'it(-[00X9A oxttsnoov 


jij/eSnxe 'X^xeuea 


Ho ‘earnsaedurei 


Figure 9>- Ambient atmospheric conditions from time of third stage ignition until fifth stage 
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Reynolds number per foot 



(b) Mach number and Reynolds number. 


Figure 10.- Concluded 
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Figure 11.- Time histories of normal and transverse accelerations 
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(a) Stations on leading edge. 


Figure 13-- Typical time histories of measured temperatures 





(b) Stations on surface center line. 


Figure 13.- Continued 
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(c) Nose apex. 


Figure 13- - Concluded 








Figure 15 . - Typical computed heating rates. 
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Distance along leading edge. In, 

Figure 16.- Leading-edge heat-transfer data. 
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Local Stanton Number x iq5 



(c) Method 1; t = 89. 0 to 89.3 seconds. 
Figure 17 •- Continued. 
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(d) Method 2 ; t = 87 . 0 to 88.6 seconds 
Figure 17.- Continued. 



Local Stanton Number x 10? 


NACA RM L58G03 



• • • • 


• • • • 

C • • • • • 


47 


Thermocouple locations 



Flagged symbols are measurements on 
other faces of pyramidal nose 
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(e) Method 2; t = 89. 0 to 89.3 seconds. 


Figure 17 . - Concluded 
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Figure l8.- Concluded. 
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